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Abstract 

In this note the Polyakov equation [Phys. Rev. E 52 (1995) 6183] 
for the velocity-difference PDF, with the exciting force correlation 
function n(y) ~ 1 — y a is analyzed. Several solvable cases are consid- 
ered, which are in a good agreement with available numerical results. 
Then it is shown how the method developed by A. Polyakov can be 
applied to turbulence with short-scale-correlated forces, a situation 
considered in models of self-organized criticality. 
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1. In Ref. [1] a method for calculating the velocity-difference PDF for 
Burgers' turbulence, based on the conjectures on OPE, Galilean and scaling 
invariance, was proposed. Starting with the Burgers' equation with a random 
Gaussian stirring force 

u t + uu x = vu xx + f(x,t) , (1) 
(f(x, t)f(x', t')) = k(x - x')5(t - t') , (f(x, t))=0, 

it was obtained that the characteristic function for the velocity-difference 
PDF (Z-f unction), determined as 

Z(n, y) = (expfaHx + y/2) - u(x - y/2)})) , (2) 
obeys the following differential master equation: 



Here n(y) is the correlation function of the external force. The \x dependence 
of the a anomaly must be chosen to conform to the scaling invariance and 
can be different depending on the scaling properties of the force correlation 
function. If for a large-scale-correlated force this function can be expanded 
as K,{y) ~ 1 — y a , then the a anomaly must depend on ji as follows: a(fj) = 
afi a , a = Using the scaling ansatz Z(/j,,y) = $(/iy 7 ) , 7 = one 
can rewrite the equation (3) in the form 



7X$" + 7(1 - 26)$' - x 2 $ = ax a § , (4) 

where x = /ly 7 . 

The unknown parameters a and b should be determined from the main 
requirement that the PDF be a positive, finite and normalizable function. 
Other possible restrictions for the theory are discussed later. 

In this note we first analyze equation (4) for arbitrary a in the force 
correlation function and consider the following special solvable cases: a — 2, 
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a — 1, a — 1/2. We then consider an application of the method of [1] to 
special kinds of external forces - forces correlated at short scales. The latter 
situation is studied in various models of self-organized criticality and has 
wide applications [2]. 

2. To begin with, we write down the asymptotics of the solutions of (4) 
for small x: 

~ H a ^— x^ + cx 2b + ... , (5) 

1 — 267 

and for large positive x: 
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<J>(x) oc exp x 2 , (6) 

where a, b, and c should be determined from the conditions mentioned above. 

We note that the most restrictive condition, the condition of normaliz- 
ability of the PDF, can be reformulated directly in terms of the ^-function. 
Indeed, the function $ must be analytical in the right half of the complex 
plane Rex > 0, and must vanish for x — > p±ioo, p > 0. This, along with the 
condition of normalization $(0) = 1, gives the quantization rule for a and b. 

Let us denote the Fourier transform of $(#) as w(z), the velocity distri- 
bution function being w(u,y) = wiu/y 1 ) /y 1 '. The integral representation for 
w(z) is: 



p+ioo 

w(z) = J e- xz ${x) dx . (7) 

p—ioo 

The asymptotics of w(z) for large positive z is determined by large x 
and is given by w(z) oc exp(— 7^ 3 /3). To find the asymptotics for large 
negative z, we deform the tails of the integration contour to coincide with 
the negative real axis. Since e~ xz decays rapidly as x — > —00, the asymptotics 
is determined by the leading singularity in the expansion (5). 

In general, two cases are possible. If < 2b, the asymptotics is 

w ~ z~ l ~~ . Such behavior is observed in numerical simulations [3], which 
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indicates that this inequality usually holds, and the 6-anomaly does not af- 
fect the asymptotics. For 2b < the asymptotics should in general be 
determined by the b anomaly, w(z) ~ z^ 2 ^ 1 , if 2b is not an integer. 

The asymptotics (5) also shows that there exist two degenerate cases. 
These are the cases when is an integer, and the corresponding term does 
not contribute to the integral. These cases (a = 1/2, 2) are solvable and 
will be considered below. We will also consider another solvable case, with 
a — 1. 

3. We start with the case a — 2. Let us Laplace-transform Eq. (4) to get 
an equation for the probability distribution w(u,y) = w(z)/y : 



w" + z 2 w' + (1 + 2b)zw = -aw , (8) 

where z = u/y, assuming the notation of [1]. All derivatives in this equation 
are with respect to z. Below we consider only the function w and drop the 
tilde sign. Asymptotics of the solution at \z\ — > oo can be easily found from 
Eq. (8): 



w oc e 2 /3 , w ~ -£-r . (9) 
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z 
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We are looking for a physically reasonable solution, with the asymptotics 



woce z3 / 3 , z — > +oo , 

W ~ z -»■ -oo . (10) 

For the w-function to be normalizable we should consider only b > 0. 
Upon writing w = \l/e~ 2 / 6 , we exclude the first derivative from Eq. (8) and 
get the Schrodinger equation for the ^-function, 



mentioned in [1]. The ground state of this equation is a positive and normal- 
izable function. This is the only solution satisfying the general requirements 
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for the PDF. Thus, for any 6 > we find the PDF as the ground state of the 
potential (11), a being the energy of the ground state. Note that the case 
6=1/2 corresponds to the left tail of the PDF ~ 1/w 2 , and the PDF obtained 
as a solution of (11), fits well the numerical observations [3] (see Figures). A 
numerical estimate in this case gives for the a anomaly a ~ 0.354. 

An important remark should be made here. Integrating Eq. (8) from — oo 
to +oo for the case 6 > 1/2, we get: 

(26 -I) J zw{z) =-aJ w{z) . 

We would like to stress that this expression does not contradict the re- 
quirement / w(u,y)udu = 0. A significant contribution to the latter integral 
can come from nonuniversal tails of the distribution function, not described 
by Eq. (8). These tails are due to spontaneous breakdown of the Galilean 
symmetry [1]. This fact should be taken into account when one compares 
the theoretical results with experimental observations. 

Nevertheless exists for which (z) =0, which corresponds to a = 0. 

To consider it, we set a = in (8), and by the substitution s = —z 3 /3 arrive 
at the degenerate hypergeometric equation: 



sw" + (7 - s)w' - aw = , (12) 

with parameters 7 = |, a = §(26+ 1).[] The positive, finite and normalizable 
solution for this case has been found in [1] . This solution can be constructed 
in the following way: the only solution, exponentially decaying at s — ► — 00 
and having power-like asymptotics at s — * +00, has the form 



w (s) = J e* (t - s)~ a t a ~l dt , s < , 
—00 

(0+) 

w (s) = J e * (t - s)- a t a ~3 dt , s > , (13) 

—00 

where in each integral the contour of the integration starts from —00, goes 
around only one of the two singular points (denoted as the upper limits) in 

1 Please, do not confuse these parameters, used only in the analysis of (12), with the 
parameters a and 7, introduced in (4). 
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a positive direction and ends up at — oo again. One of these solutions can 
be analytically continued to the other one only if a = n — 1/6, where n is 
any integer number. It is interesting to note that this exact quantization 
rule can be easily obtained as the Bohr-Sommerfeld condition for quantum 
mechanics considered above, with zero energy.^ Positivity of the solution 
requires n — 1. 

For the other degenerate case, the correlator of the external force has the 
form: n(y) = 1 — y 1 ^ 2 . This force leads to a differential equation for the 
w(z)-function, analogous to Eq. (8): 

w" + -z 2 w' + (- + b)zw = aw' , (14) 

where z = u/y 1 ^ 2 . Asymptotics of the left tail of the solution is given by (10). 
Excluding the first derivative from this equation, we obtain the Schrodinger 
equation for the function \1/ = w exp(z 3 /12 — az/2): 

_ + ( z ! _ _ bz ) $ = ^ . (15) 
V 16 4 J 4 v ' 

As in the previous case, one can find the solutions as the ground states 
of this equation. The numerically observed PDF [3] has the left tail ~ 1/u 3 
in the considered case. The same PDF can be obtained from our equation if 
we set b = 1, i.e. when the (3 anomaly is absent. One can then numerically 
obtain a ~ —0.473. A comparison of the whole PDF with the numerical 
results [3] shows a very good agreement (see Figures). 

Equation (15) has another interesting feature. For sufficiently large a the 
potential is a two-well function, and one can show that the ground state is 
concentrated in the right well. Rescaling the variable z — > zar l l A , one can 
expand the potential near the bottom of this well to obtain: 

- + ^z 2 ^ = bV2^ . (16) 

It is interesting that a has dropped out of this equation. Energy of the 
ground state is l/v2, which gives 6=1/2. A numerical estimate shows that 

2 This was pointed out by V. Gurarie 
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this result holds accurately already for a > 1.5. One could say that with 
decreasing energy, the system freezes at the point 6 = 1/2. For this value 
of 6 Eq. (14) can be solved exactly. 

To analyze the last case, a = 1, let us work in the x-representation. Note 
that by the substitution £ = x 3//2 one can cast equation (4) into the form: 

|c^+(|-26)^-(C + a)$ = 0. (17) 

This equation can be solved by the Laplace transform. The solution with 
the correct asymptotics is: 



(0+) / r- \oa 

$( x ) = Cx-H ai+a2+ \Vi* 3/2 I e T r Ql r + 2W|x 3/2 dr . (18) 

with 

«i = -- [1 + 46/3] + a/VQ , a 2 = -- [1 + 46/3] - a/VQ ■ 

<&(x) will be an analytical function for Rex > 0, and a decaying function 
for x — > p±ioc only when ai = n or «2 = m, where n is any negative integer 
number and m is any non-negative integer number. The only possibility of 
getting $(0) = 1 is ai = n, which gives the following quantization rule: 
a/v^6 — 26/3 = n + 1/2. Positivity of the solution forces us to select n — — 1, 
and (18) reduces to $ = exp y^|a;l .0 

4. We now show how, using the ideas of [1], one can consider turbulence, 
excited by a force correlated at small distances. We assume: 



(f(x, t)f{x>, 0) = D(L, r c )S{t - if) , x = x' , 

(f(x,t)f(x>,t'))=0 , |x-x'|»r c , (19) 

where r c is the correlation length, L is the dimension of the system. We are 
interested in the velocity-difference PDF for large distances: y ^> r c . We shall 

3 V. Gurarie informed me that this solution can probably be obtained by the instanton 
method, along the lines of [4]. 
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assume that the integral of the force correlation function / (f(x, t)f(x', t'))dx 
remains finite as r c — > 0. In previous treatments of this problem, based 
mainly on the renormalization group approach, the limit r c — > was taken 
before considering small v [2]. In fact, v was kept finite there and entered 
the final answer. Here we first take the limit v — > 0, and then consider the 
second limit r c /y — > 0. The answer should not depend on either v or r c . 

Assuming the Galilean and scaling invariance we obtain the following 
equation for the Z function: 



Now one can take the limit r c — > 0. For this purpose consider the singularity 
in the coefficient D: 



D(r e ,L) = ±- + 0{l), r c ^0.(^y . (21) 

To get a finite limit, the singularity in the anomaly a should cancel the 
same singularity in D. Assuming such a cancellation we will use the same 
letters D and a for the nonsingular parts of these coefficients. Introducing a 
new function <3>: Z(fx,y) = ^(fiy 1 ^ 3 ) we obtain an equation for the Laplace 
transform of <£>, completely analogous to the Eq. (8): 



3(D + a)w" + z 2 w' + (1 + 2b)zw = , (22) 

where z = u/y 1 ^ 3 . Note that except for the cancellation of the singularity 
in the D- coefficient, the a-anomaly does not play an important role here. 
Introducing the new variable s = — 9 p +a ) , we get exactly equation (12). This 
equation has a unique solution, which corresponds to the velocity-difference 
PDF with the left tail l/u§. 

Let us discuss the conditions for the Galilean and scaling invariance as- 
sumed in the approach. The force (19) generates finite v rms : 

v rms = lim . [r c D(L,r c )] 1/3 , 
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and the Galilean invariance holds for u <C v rms . The condition for the scaling 
invariance is y <C L. 

5. Finally, we discuss an important general restriction that can be im- 
posed on the theory. This follows from the physical condition of positivity 
of dissipation and was proposed by A. Polyakov.0 It can be obtained if one 
notes that the operator 



rj2 



dx 2 

is not singular if x does not coincide with any other Xi. Therefore 



r%2 

lim v— e Xu{x) = , (24) 
u^o dx 2 



which leads to 



a(X)Z + ^Z = -hmu(X\le^ + ---) , (25) 

where we use the notation of [1]. The r.h.s. of this expression is negative. 
The function a (A) is analytical in the right half of the complex plane, and 
may have a discontinuity at the imaginary axis. Summing up corresponding 
expressions for Ai = /x/2, X\ = x + y/2 and A2 = — A*/2, xi = x — y/2, we get 
the following necessary condition, that must be valid for all non-negative x: 



ax a $ - 2(1 - &)$' < , (26) 

where b — $ + 1. 

One can easily see that this condition is rather strong and allows one to 
considerably restrict the possible solutions of (3). For example, it prohibits 
the solutions with b < 3/4 for the case a = 2 and, probably, forces the (3 
anomaly to vanish for a < 1/2. 

Nevertheless, one can see that this inequality is absent (or, at least, the 
above arguments do not work) if we consider the dissipation in the form 

4 A. Polyakov, 1996, unpublished. 



9 



ud 2p /dx 2p , with p > 1 (the so-called hyperdissipation). This is the case for 
which the numerical simulations [3] have been performed. One can expect 
that the structure of the shock front will be changed in this case. For example, 
one can show that for the dissipation of the form v\u xx — Vz u xxxx there exists 
the stationary solution, proportional to the linear combination of — tanh(x) 
and — tanh(x)/ cosh 2 (a;). 

The intriguing question whether these different dissipations lead to dif- 
ferent stationary regimes of the Burgers' turbulence and to what extent the 
found solutions can correspond to them has yet to be answered. 
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Figure captions. 



The following results of numerical simulations, Figures 3 and 5, are taken 
from the paper by V. Yakhot and A. Chekhlov [3]. Fig. 3 shows the collapse 
of the PDFs in the universal region of Au, and in our notation corresponds 
to a = 1/2 in the force correlator. Fig. 5 shows the collapse of the PDFs 
in the universal region for a — 2. For comparison, we have depicted the 
theoretical results by dots. 



11 



